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We develop a theoretical proposal for the charge-Kondo effect in mesoscopic normal- 
superconductor-normal heterostructures, where the superconducting gap exceeds the electrostatic 
charging energy. Charge-Kondo correlations in these devices alter the conventional temperature- 
dependence of Andreev reflection and electron cotunneling. We predict typical Kondo temperatures 
of > lOmK, and suggest experimental architectures that combine superconducting charge-qubits 
with semiconducting nanowires at cryogenic temperatures. 



PACS numbers: 



I. INTRODUCTION 



Five decades after its discovery, the Kondo effect 
is regarded as an archetype many-body phenomenon 
that interconnects strongly correlated systems such as 
atomic nuclei, heavy fermion compounds, superconduct- 
ing cuprates and quantum dot devices ^ The Kondo ef- 
fect emerges from the interaction of localized, degenerate 
states with itinerant degrees of freedom. It is character- 
ized by a low-temperature infrared divergence in pertur- 
bative calculations of physical observables such as resis- 
tivity and magnetic susceptibility. Often the localized 
degrees of freedom originate from electronic spins and 
their degeneracy is tied to the usual spin degeneracy. 
The ensuing spin-Kondo effect (SKE) is realized when di- 
lute magnetic impurities arc embedded in a non-magnetic 
host, or when conducting electrodes are coupled to quan- 
tum dots<2 

Aside from spin, any degenerate two-level system that 
is coupled to a fermionic bath is a potential host for 
the Kondo effect. In particular the charge counter- 
part of SKE, known as charge-Kondo effect (CKE), 
arises when a Fermi sea hybridizes with dilute impurities 
containing two degenerate charge states. This variant 
has received sustained theoretical attention for the last 
two decades, being investigated in non-superconducting 
single-electron devices^— as well as in valence-skipping 
elements with attractive onsite interactions ("negative- 
U molecules" ) .— ~— In the former case the electrostatic 
energy can be engineered by gate electrodes to be the 
same for two states that differ by one electron, whereas 
in the latter case chemistry dictates that the charging 
energy be degenerate between states that differ by two 
electrons. Incidentally, CKE in single-electron devices 
constitutes a paradigm of the two-channel Kondo effect, 
where the channels originate from spin. 

As opposed to SKE, which has been thoroughly 
observed both in its single-channel and two-channel 
versions,— the experimental detection of CKE remains 
at a primitive stage. Experimental challenges abound 
for the implementation of CKE in non-superconducting 
devices. On one hand, CKE is sensitive to and washed 
out by background charge noise, much like SKE is sen- 



sitive to and washed out by magnetic fields. On the 
other hand, both the charge-Kondo temperature and 
the charging energy of the system arc required to be 
large compared to the single-particle energy-level spac- 
ing in the grain, which is difficult to achieve in semi- 
conductors. Prospects may be better for metallic sys- 
tems, wherein single-particle energy levels form a near- 
continuum. Yet metallic devices often involve junc- 
tions with multiple channels, and unfortunately the 
charge-Kondo temperature in non-superconducting sys- 
tems scales exponentially unfavourably with the num- 
ber of channels.— Some of these problems might be cir- 
cumvented by using atomic point contacts^ or else by 
resorting to resonant tunneling devices i 16 ' 17 Altogether, 
attempts to measure fingerprints of the CKE in non- 
superconducting systems have thus far met with sugges- 
tive yet inconclusive outcomeSi 18 ' 19 

Recent experimental efforts 2 ^ concerning negative-U 
molecules appear to have met with more success: charge- 
Kondo temperatures of ~ 5K have been reported in PbTc 
doped with valence-skipping Tl. Unfortunately, real ma- 
terials with valence-skipping compounds are relatively 
rare, more complicated and less tunable than artificial 
single-electron devices. 

In this work we theoretically demonstrate that CKE 
can also occur in artificially fabricated mesoscopic 
superconducting islands that are connected to non- 
superconducting leads (NSN devices). We concentrate 
on superconducting grains whose many-particle energy- 
gap (A) exceeds the electrostatic charging energy (E c ). 
These grains behave somewhat like giant negative-U 
molecules, and they are to valence-skipping elements 
what quantum dots are to local magnetic moments: more 
easily manufacturablc and more readily controllable sys- 
tems, albeit at the price of a paramctrically lower Kondo 
temperature. 

Our proposal is a natural and perhaps obvious exten- 
sion of Refs. JjIJ] and [l3|, as it combines attractive 
(phonon mediated) and repulsive (electrostatic) Coulomb 
interactions with a dense single-particle energy spectrum. 
Surprisingly, there have been no thorough studies of CKE 
in superconductors with A > E c . The relevant literature 
is limited to some peripheral statements ; 21 ' 22 along with 



a tacit assumption that CKE in superconducting dots 
is conceptually similar to CKE in non-superconducting 
single-electron devices (NNN devices). Contrary to this 
view, our results aim to draw the attention of theorists 
and experimentalists alike towards the study and search 
of charge-Kondo correlations in NSN systems. 

The rest of this paper is organized as follows. In Sec- 
tion [il] we introduce the basic microscopic model for 
a NSN heterostructure. The island is assumed to be 
chaotic and smaller than the superconducting coherence 
length, thereby allowing for coherence between electrons 
tunneling through different junctions. In Section IIIII we 
map the superconducting island onto an artificial spin 
1/2 describing two charge states. This mapping is well- 
known and has already been exploited in existing Cooper 
pair boxes. In Section IIVI we attach normal metallic 
leads to the superconducting island, and find that at low 
enough energies and near the charge degeneracy point 
this system is described by an anisotropic Kondo model. 
The degree of anisotropy and the magnitude of the Kondo 
couplings can be changed by tuning the ratio A/E c . 
In Section [V] we provide quantitative estimates for the 
charge-Kondo temperature in experimentally realizable 
NSN devices. In addition, we present a detailed discus- 
sion on how our proposal for CKE fares in comparison 
with previous proposals concerning negative-U molecules 
and non-superconducting devices. In Section [VII we cal- 
culate the fingerprints of CKE in physical observables of 
NSN systems. Specifically, we focus on the temperature- 
dependence of the zero-bias conductance at low tempera- 
tures. In Section lVlIl wc determine the fate of CKE when 
the normal-superconducting junctions arc highly trans- 
parent; the outcome depends on the strength of electron- 
electron interactions in the normal metallic leads. Sec- 
tion [VlTT] is devoted to a short summary and conclusions, 
and the Appendix contains technical calculations that 
rigorously justify the considerations of Section IIVI 



II. MICROSCOPIC MODEL 

Let us suppose we have a mesoscopic superconducting 
grain that is weakly coupled to normal metallic leads. Its 
Hamiltonian can be expressed as 



U = Ui + Ha+Ut 



^ C cxka Cako 



Ut 



aha 

" endLd™ + E C (N - N g f + T0f 



/ j ^n^ na ^na 



/ , ta kn C akcr d n( j 
akan 



h.c. 



(1) 



%\ is the Hamiltonian of the leads. c' aka creates an elec- 
tron with momentum k and spin er, in a lead labeled by 
a (a = L(cft), R(ight)). £& is the energy dispersion of the 
itinerant fcrmions measured from the Fermi energy. For 



now we assume that there is only one conduction chan- 
nel on each lead; as such, k is the momentum perpen- 
dicular to the interface between the normal-metal and 
the superconductor (NS interface). Generalizations to 
multi-channel leads, which are necessary when the linear 
dimensions of the NS interface exceed the Fermi wave- 
length, will be discussed in Section V. Likewise we ne- 
glect electron-electron interactions in the leads, although 
this assumption will be relaxed in Section VIII. 

Hd is the Hamiltonian of the superconducting quan- 
tum dot. S na creates an electron in the n-th single- 
particle level of the island. The corresponding energy 
(measured with respect to the Fermi energy) can be writ- 
ten as e„ = nS, where 5 is the average single-particle level 
spacing in the dot. Hence, in this model every energy- 
level is non-degenerate, save for the twofold spin degen- 
eracy. N = J2na d nadna and f = Y^ n rf n-lAt are the 
number and pairing operators in the island, respectively. 
7] is the effective (phonon mediated) coupling between 
electrons in the island. Hereafter we adopt the BCS 
approximation^ by introducing an s-wave superconduct- 
ing gap A = rj(T). The BCS approximation is accurate^ 
provided that A > S, which constitutes the regime of in- 
terest in this paper. For reasons exposed below, we take 
A > E c , where E c = e 2 /2Cs is the charging energy of the 
island. Cs = Cl + Cr + C g is the total capacitance de- 
vice, Cl(b) is the capacitance of the left (right) NS junc- 
tion and C g is the gate capacitance. In mesoscopic NSN 
heterostructures most of the capacitance emanates from 
the junctions and scales roughly as Cs ~ eA/d, where e 
is the dielectric constant of a thin insulating layer that 
separates the normal metal from the superconductor, A 
is the surface-area of each NS junction and d ~ lnm is 
the thickness of the insulating layer. N g is the charge 
induced in the island by the gate electrode. We neglect 
intra-dot exchange interactions, which is adequate given 
the spin-singlet order parameter of the superconductor. 

7-Lt models the single-particle tunneling between the 
island and the leads, t akn being the tunneling amplitude. 
Throughout this paper (with the exception of Section 
VIII) we treat the tunneling term as a weak perturbation. 
Hence we require that 5 » T, where T is the broadening 
of the energy-levels of the island due to tunneling between 
the island and the leads. 

Our model relies on precepts from random matrix the- 
ory (RMT) and is valid only for dots without spatial 
symmetries and for states in the vicinity of the Fermi 
energy ££5- According to RMT, single-particle states with 
|e„| < Et are non-degenerate except for their spin degen- 
eracy, where Et is the Thoulcss energy. Et = hD/L 2 
for a dirty sample (where D is the diffusion constant 
and L is the linear dimension of the superconductor) and 
Et = hvF I ' L for a clean sample (where vp is the Fermi 
velocity). Anticipating the fact that charge-Kondo cor- 
relations originate mainly from states with \e n \ < A, we 
request Et > A so that Eq. (fl]) can capture CKE quan- 
titatively. This condition is tantamount to L < £, where 
£ is the superconducting coherence length. In sum, the 



hierarchy of relevant energy scales reads 
E T > A > E c > 5 > T. 



(2) 



L < C implies electron coherence between the two junc- 
tions, which in turn gives rise to elastic electron cotun- 
neling and crossed Andreev reflection (Section IVI[) . Un- 
mentioned as it is in Eq. @, we note that the Dcbye 
frequency u-q (such that A ^ iff |e„| < ujjj) is much 
larger than Et for all cases discussed in this paper. 



III. ISOLATED SUPERCONDUCTING ISLAND: 
A PSEUDOSPIN 1/2 QUANTUM IMPURITY 

In this section we ignore the leads and simply con- 
sider an electrostatically gated superconducting island, 
which contains a definite number of particles TV because 
[TV 2 , TV] = [ftf, TV] = 0. Its ground state depends on TV 
as2£ 

E (N) = E C (N -No- TV 3 ) 2 + |(1 - (-1)"), (3) 

where TVo is the number of electrons in the island when 
the gate voltage is zero. We assume without loss of gen- 
erality that TVo is even so that TV — TVo and TV have the 
same parity; herein we absorb TVo hi the definition of TV. 
The last term of Eq. ([3]) imposes an additional energy 
cost whenever TV is odd, because of a single unpaired 
quasiparticle that must reside outside the superconduct- 
ing condensate. The lowest energy this quasiparticle can 
have equals A. If A > E c and if the gate voltage is tuned 
to N g = 2M + 1 (for any integer M), the lowest energy 
eigenvalues are (ordered from lower to higher energy) 

E(2M) = E(2M + 2) = E C 
E(2M + l) = A 

E(2M - 1) = E(2M + 3) = AE C + A 
E{2M + 4) = E(2M - 2) = 9E C 

(4) 

where the ordering of the states changes if E c < A/5. 
The fact that E(2M + 1) < (E(2M - 1), E(2M + 3)) will 
be important for establishing CKE in Section HV1 

Since the ground state of the island is doubly degen- 
erate, it behaves as a spin 1/2 system at low enough 
temperatures (T < T*). For later purposes we introduce 
a pseudospin label: 



\2M) = \U) 5 \2M + 2) = \H ), 



where the subscript simply labels the location of the 
island. From now on pseudospin and real spin degrees 
of freedom will be denoted by double and single arrows, 
respectively. This qubit is well-separated from a dense 
forest of many-body states, which have excitation ener- 
gies greater than A — E c . Due to entropic issues 2 ^ that 



arise from having an approximately manifold degenerate 
set of excited states, 



rp ■■■-■ 



A- E c 



(0) 



is smaller than the excitation energy gap A — E c . 

The mapping to a two-level system holds even when 
M is macroscopically large, and has been amply corrobo- 
rated by experiments<2Zr— The degeneracy between \2M) 
and \2M + 2) may be lifted by tuning the gate voltage 
away from the degeneracy point, which is akin to apply- 
ing a pseudospin magnetic field along z, or by coupling 
the island to a superconducting reservoir, which acts as 
a pseudo-magnetic field in the xy plane. The charge- 
Kondo effect discussed below emerges when these pseu- 
dospin magnetic fields are very weak. 

For completeness, we note that there is a twofold 
charge-degeneracy point even when A < E c . In this case 
the degeneracy occurs between states that differ by a sin- 
gle electron^ and therefore this problem is adiabatically 
connected to that of Ref. [3| . 



IV. COUPLING TO THE LEADS: 
MODEL 



KONDO 



In this section we consider the influence of conducting 
leads on the superconducting island with charge degener- 
acy, at temperature T < T* . We assume all non-thermal 
perturbations to be small compared to A — E c Then the 
NSN heterostructure behaves as a localized pseudospin 
1/2 interacting with a bath of itinerant fermions, and 
Eq. (fTJ) can be truncated on symmetry grounds into an 
effective Hamiltonian 



+ E[ A r' 52s -'+ A " Q ' (s+§- a , +h.c. 



,(7) 



where 2S° = | fto) (ft, 1 + 1 h) (h \,2S* = \ to) (fto |-| ^o 
) (JJ-o | and S + = \ ffo) (-IJ-o I are the pseudospin operators 
characterizing the state of the superconducting island, 
and 



E 






Ca' k'' 



(8) 



kk f 



denotes the pseudospin density in the leads (i € 
{0, x,y, z}). In addition, t is a vector of Pauli matri- 
(5) ces (r is the identity matrix) and c a k a is a pseudospin 



operator defined as 



-afct 



-afcf i C-ak^ 



-J 



"a,-fe4-' 



(9) 



Physically, s is the spin density along z, s z is the charge 
density and s is the pairing operator. In Eq. ([7]). the 
term involving S° (S z ) describes normal spin-dependent 



(spin-independent) scattering and the term oc 5 ± de- 
scribes Andreev processes whereby the number of Cooper 
pairs in the island changes by one; the xy symmetry 
is simply a reflection of gauge invariance. Ignoring 
real-spin- flips (which are unimportant at low energies), 
Eq. (|7J) is the most general form of exchange interac- 
tion between the localized and itinerant pseudospins. We 
have omitted perturbations that break charge degener- 
acy; these can be easily incorporated in Eq. ([7]) as pseu- 
dospin Zeeman fields. 

Am ^ are the Kondo couplings. Their microscopic ex- 
pressions can be extracted perturbativcly from Eq. ([T]) 



(/|«-ff|») 



v (f\H T \n)(n\H T \i) 
^ E, - E n 



(10) 



where |i) (|/)) is the initial (final) state with energy Ei 
(Ef) and \n) is an intermediate state with energy E n (see 
Table IJ). It follows that 

, QQ ' y^ (it ak ;ito \n T \n)(n\H T \ iWilro) 

A o - 2^ 



\ aa 

A n 



E 

n 

E 



E 1tck;ito ~ E » 
(1\ a kAo\UT\n}(n\H T \1\ a >k>Ao) 

E 1tch;ii-o ~ E n 

i'frak-Ao\HT\n){n\HT\'fra>k>-Ao) 



E. 



iW;lk> 



E 



II 



Y (1\ a k-Ao \n T \n)(n\n T \ 1\a>k>Ao) 

„ E 1t«kM 



En 



X°[ a 



Y (fak-Ao \n T \n)(n\H T \ K'k'-Ao) (u) 
„ E 1tc,k;U - E n 

where | f[- afc ) = 4, fct |0) and | ty ak ) = c QfcJ j0) represent 
spin-up electrons (spin-down holes) at the Fermi surface 
(|0) stands for the filled Fermi sea in the leads). Table U 
collects all the intermediate states \n), along with their 
energies E n . 

We proceed with the evaluation of the tunneling matrix 
elements using the BCS formalism. For instance, the 
matrix elements for an electron or a hole tunneling from 
lead a to the n-th level in the island are 

(1\ a k;1\o \H T it t \ 1\o) - (1\ a k-Ao \H T ~fl f \U)^t akn u n 

(Kk-Ao \n T il t \ 1\ ) =s Aak-A \n T ii t \U)^-ta kn vn 

where Ana = u nd na + <rv n d n ^ a is the BCS quasiparti- 
cle creation operator, u„ = (1 + e n / ' \]t n + A 2 )/2 and 
v n = 1 — w 2 are the standard BCS coherence factors (we 
assume u n and v n to be real without loss of generality be- 
cause we consider only one superconductor). Other ma- 
trix elements may be computed similarly. Even though 
intermediate states like j n(J \2M + 2) do not have a well- 
defined number of particles, they are excellent approxi- 
mations to the true eigenstates of the superconducting 
island when A >> 8. 



The explicit expressions for the couplings can now be 
read out from Table HI 
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A || 
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-3E C - y/e 2 n + A 2 E c - y/el + A 2 



\aa _ 9 \ +* +* 



V anka ' k ' n E c -^7iT^' 



(12) 



It is easy to verify that Ao = by invoking particle-hole 
symmetry for the energy spectrum of the superconduct- 
ing island. Since Am is associated with normal scattering, 
it involves virtual transitions that do not flip pseudospins. 
The sign of A|| will be important in the considerations 
below. Aj^ is associated to Andreev scattering, hence it 
involves virtual transitions that flip the pseudospins of 
both the island and the leads. The sign of Aj^ is unim- 
portant as it may be reversed by a gauge transformation. 
In principle, evaluating the sums in Eq. (fT2"|) requires 
a detailed knowledge of the tunneling amplitudes t ankl 
which in turn depend on details of the wavefunctions 
at the NS interfaces. In practice, simplifying assump- 
tions may be quantitatively satisfactory. For instance, 
t an k — tan is a good approximation for NS junctions 
that are atomically thin (recall that k is the momen- 
tum perpendicular to the NS interface). This renders 
Am^_ independent of k and k' . Similarly, we conceive 



t* t i 



tAa' and S„_t„/_ 



tit*, with the understand- 



ing that the n— dependence is weak for states that are 
contained within a very narrow strip around the Fermi 
energy. This approximation is motivated by the fact that 
the main contributions to An and Aj^ in Eq. (|12p originate 
from states with |e„| < E c < Et and |e n | < A < Ex, 
respectively. 

Under the above assumptions, the tunneling ampli- 
tudes may be taken to be real and the Kondo couplings 
read 



a; 



= A 



LL 



1 t R /t L 

t R /t L {t R /t L f 



(13) 



Eq. (|13p may be diagonalized with a unitary 



transformation^ which converts Eq. ([7]) into 
single-channel Kondo model: 

Hsck = Ux + [\\\S z st + Ajl (S+Sr + h.c.)] , (14) 
where s\ = J2 k ,k' Ak a T la'C\k' a' and 



Clfeo 

A,, 



{t L CLka + t R C Rka )Nt\ + t\ 



— \ LL _i_ \RR 

- A ||,_L + A ||,_L 



(15) 



Thus only one out of the two conduction channels is 
coupled to the superconductor. Starting from Eq. (|12[) 
and using e„ = n8 we derive explicit expressions for the 
Kondo couplings at second order in tunneling amplitude: 



v\±_ = -?f±(x) ; v\\ 



r 



/||(a0, 



(16) 



TABLE I: Virtual elastic processes to second order in single-particle tunneling (the extension to fourth order processes is 
deferred to the Appendix). The initial and final states are in the truncated, low-energy Hilbert space whereas the intermediate 
states trespass into the high-energy sector. We ignore inelastic processes because they are suppressed at T <C T* . The energy 
of all initial and final states is ~ E c because the most important tunneling events involve electrons/holes at the Fermi surface 
(&~&,~0). 
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where x = E c /A, 
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, (17) 



i/ is the density of states of the leads at the Fermi energy 
and r = ~Kv(t\ + t R ) denotes the broadening of the en- 
ergy levels in the dot due to its coupling to the leads. Y 
is related to a dimensionless parameter representing the 
transparency of the tunnel junctions: 



9t = 9l + gR = 



(t| + t R )i 



= -x « lj 



(18) 



where gi, and gn are the conductances of the left and 
right tunnel junctions in units of 2ire 2 /h. In the deriva- 
tion of Eq. (|17l) we approximated J2 n ^Y I- ^ n > * ms 
is justified in mesoscopic superconducting islands, where 
S < 10~ 3 A. In Fig. ((T|) we plot f± and /||. An essentially 
identical expression for Aj_ was found by Hekking et al. 
in Ref. [3J]; however, these authors did not contemplate 
the possibility of a Kondo effect and did not evaluate 
Am, whose positive sign is important in order to produce 
CKE. 

As demonstrated in the Appendix, Kondo correlations 
become apparent only when we carry out perturbation 
theory to fourth order in tunneling. Indeed at fourth 
order both coupling constants Aj_ and Am develop in- 
frared singularities, in a manner that is consistent with 
the renormalization group equations of the anisotropic 
Kondo model)2£i2& 



dA|| 

~df 
dX ± 

~oY 



= 2uX\ -2v 2 X\\X\ 



^A ± A||-^ 2 A ± (Af | +Ai), 



(19) 



where I = ln(T*/T) and by definition A||,_l(Z = 0) = 
Aim (T*) are the values of the Kondo couplings at an 
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FIG. 1: Parameters f±, /|| and g± as a function of E c /A. 
These parameters (defined in the main text and in the Ap- 
pendix) determine the values of the Kondo couplings at an 
energy scale T* (see Eqs. ® and pO])). When E c is suffi- 
ciently close to A (precisely how close is contingent on the 
value of r/<5) the bare couplings become large, which inval- 
idates the perturbative renormalization group approach em- 
ployed in this work. For typical mesoscopic superconducting 
islands E c /A ~ O(l) and therefore X n (T*)/X±(T*) ~ O(l). 



energy scale T*, below which the Kondo model becomes 
applicable. The present analysis is quantitatively reliable 
provided that Aim (T*) -C 1. We show in the Appendix 



that 
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where g± and g\\ are dimcnsionles functions of A/E c (the 
former is plotted in Fig. ([1])). Eq. ([2"U|) is valid provided 



that the 0(T 2 ) terms are smaller than the 0(T) terms; 
for T/5 -> Eq. flU) reduces to Eq. jTB]). The loga- 
rithmic terms in Eq. ([20]) suggest that the Kondo cou- 
plings begin to renormalize from energy scales of order 
A. Aii(T*) > ensures that both An and A^ will flow 
to strong coupling through Eq. (fT9|) . The energy scale at 
which max{i/Aj_(Z c ), i/Am(7 c )} ~ 1 is given by the Kondo 
temperature: 



T K = T*exp(-Z c ), 



(21) 



where T* should be replaced by 8E C when E c < T* 
(recall Eq. (gj). 
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V. DISCUSSION 
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In this section we derive quantitative estimates for the 
charge-Kondo temperature in NSN heterostructures, ex- 
amine the influence of multi-channel junctions, propose 
various experimental NSN platforms wherein CKE might 
be observed, and compare our proposal of CKE in NSN 
heterostructures with previous proposals of CKE in other 
systems. 



A. Estimates for the Charge-Kondo Temperature 



FIG. 2: Kondo temperature (Eq. ([21}) (in units of A) as a 
function of E c /A, for T/5 ~ 0.08. We choose A/S ~ 1000, al- 
though the results are only weakly sensitive to A/S. For this 
set of parameters, the Kondo temperature can reach ~ lOmK 
by selecting E c /A in the appropriate range. The decrease 
in T* as E c /A increases is overcompensated by the gain in 
the bare Kondo coupling. In our estimates we have ignored 
both g± and g\\ . The perturbative analysis breaks down when 
z/A||,x(T*) ~ O(l). For T/S ~ 0.08 this breakdown is appar- 
ent at -E c /A > 0.5. Inset: Kondo temperature (in units of A) 
as a function of r/<5, for E c = 0.5A. 



The Kondo temperature of our system can be quanti- 
fied combining Eqs. ([16 ]) . ([17 ]) . ([15 ]) . ([20 ]) and J2IJ. For 
pedagogical purposes we begin by discussing a few par- 
ticular cases that can be solved analytically. 

The first special case corresponds to E c ~ 0.5A, for 
which Aj_(T*) ~ A||(T*) and 

T K (A^2E C ) ~ T*cxp(-1/2^A ± ) 

~ Acxp(. g± )exp(-^j, (22) 

where for simplicity we have neglected third order terms 
in the RG equations. In the second line of Eq. ([2"2"| 
we have made a Taylor expansion under the assumption 
that (T/S)g±/f± < 1 and 2z/Aj_ ln(A/T*) < 1. For 
mesoscopic superconductors, ln(A/T*) < 2 and hence 
the latter condition requires that T/5 < 0.1. A typical 
Kondo temperature in this case (using f± ~ 1.4) reads 
Tjc ~ 10 _3 A ~ 1 — lOmK for conventional superconduc- 
tors. It is interesting that even though the Kondo Hamil- 
tonian is valid at energy scales below T* , the prefactor of 
the Kondo temperature in Eq. (|22[) is effectively shifted 
from T* to A. An analog situation occurs between E c 
and 5 in SKE^l 

Another simple limit is that of E c <C A, in which case 
A||(T*) < X±(T*) and accordingly^ 



T K (A > E c ) ~ T* exp(-7r/4i/A_L), 



(23) 



where once again we have neglected third order terms in 
Eq. Unj). The exponent of Eq. (f2"3"]) is a factor 7r/2 more 



negative than the exponent of Eq. ([22]) , and f±(E c <C 
A) ~ 1 < f±(E c ~ 0.5A). Consequently, for a given 
value of T/5, T K (A > E c ) < T K {A ~ 2E C ). 

A third special case consists of A||(T*) >> Aj_(T*). 
However, in our model this regime does not arise for any 
value of E c /A and thus we disregard it. 

We now solve more general cases numerically. Fig. ([2]) 
displays the charge-Kondo temperature for a generic 
range of parameters. For fixed T/5, Tk increases rapidly 
with E c /A (at least up until E c /A ~ 0.6). Our findings 
establish Tk — lOmK for superconductors such as Al or 
Ta, and Tk — 20mK for superconductors such as Nb and 
Pb. 

As -E^/A — > 1, higher order tunneling processes (en- 
coded in <7|m ) become important, u\j_ n (T*) is no longer 
small and correspondingly our perturbative calculation 
of Tk becomes unreliable. The divergence of g± towards 
negative values (Fig. [T]) casts some doubt on whether the 
Kondo temperature will continue to increase as E c /A — > 
1. 

In the inset of Fig. ([2]) we illustrate the dependence 
of Tk on T/5. For a given value of E c /A, Tk increases 
exponentially with T/5. Only small values of T/5 are 
allowed in our calculation, which treats tunneling as a 
weak perturbation. In the strong tunneling regime (r 3> 
5) the charge of the superconducting island fluctuates 
strongly and it is a priori unclear what the fate of CKE 
will be. This issue will be addressed in Section IVlIl 



B. Influence of Multichannel Leads 

Thus far we have been considering only one conduction 
channel per lead. This is appropriate for atomic point 
contacts, wherein itinerant electrons are fully character- 
ized by their momentum in the direction perpendicular 
to the NS interface. Nonetheless, due to the small Fermi 
wavelength of metals, most NS contacts contain multi- 
ple channels that describe e.g. transverse momenta of 
electrons in the lead. Multichannel leads can be read- 
ily incorporated into Eq. (JlJ by rewriting the tunneling 
Hamiltonian and the Hamiltonian for the leads: 



kai 
Ht = ^ t knicl ai d n(J + h.i 



(24) 



kani 



where i = 1, ..., 2N is a channel index that includes the 
which-lcad label. In wide junctions the tunneling am- 
plitude depends on the channel index i. In contrast, we 
have taken the kinetic energy of the itinerant fcrmions to 
be independent of i; this can always be ensured by proper 
rescaling of the itinerant fermion operators. In addition, 
we have assumed that each energy level in the supercon- 
ducting dot remains only twofold degenerate, due to the 
chaotic motion of electrons in the dot. 

In principle it is possible that only one linear com- 
bination of the channels couple to the superconductor. 
This is the case when any one of the following assump- 
tions (from more to less stringent) is satisfied at energies 
|e„| < A: (i) t^ni — ifci (energy- independent magnitude 
and sign of the tunneling amplitude), (ii) tkni — drfiu 
(separable tunneling amplitude), (hi) tknitk'nj — tkitkj 
(energy-independent magnitude of the tunneling ampli- 
tude, and coherent tunneling through different junctions 
due to L < (). With this proviso one can always make a 
unitary trasformation in the spirit of Eq. (TT5|) and recover 
an effective one-channel Kondo model. 

When the assumptions (i)-(iii) above fail, a multichan- 
nel Kondo effect invariably follows. In NNN devices Tk 
drops exponentially^ when the number of channels in the 
normal leads and in the normal metallic grain is much 
larger than one. This argument does not apply to NSN 
devices insofar as the superconducting island is describ- 
able by random matrix theory and L < £. However, 
even in our case Tk is generically lower for the multi- 
channel than for the single-channel Kondo model. Sup- 
pose for instance that there arc N ^> 1 separate Kondo 
couplings of similar magnitude. In order to maintain the 
Coulomb blockade, each bare Kondo coupling must scale 
like ~ 1/-/V, i.e. each Kondo coupling must be very small. 
Eq. (TK)]) then dictates that the Kondo temperature will 
be very low. This state of affair changes whenever a uni- 
tary transformation can map multichannel leads onto an 
effective single-channel lead, because in that case the only 
nonzero Kondo coupling in the transformed basis scales 
as the sum of all the original couplings (see Eq. (|T5j) for 
a N = 2 example of this). 



In order to arrive at a single-channel Kondo problem in 
presence of wide NS junctions it is helpful that the Kondo 
couplings be of second order in the tunneling amplitude, 
because condition (iii) above is more realistic than ei- 
ther (i) or (ii). In case of CKE in non-superconducting 
dots£ the xy-component of the Kondo coupling is ~ 0(t), 
whereas the z— component is ~ 0(t ). In such situation 
one is generally left with a multichannel Kondo problem. 
This is a potentially important difference between CKE 
in NSN and NNN heterostructures. 



C. Possible Experimental Setups 
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FIG. 3: Two possible arrangements for the experimen- 
tal observation of the charge-Kondo effect in NSN sys- 
tems. S=superconducting island, G=gate, N=normal metal, 
NW=semiconducting nanowire. (a) Metallic-semiconducting 
hybrid setup (top view). The width of the semiconduct- 
ing nanowires is comparable to their Fermi wavelength, 
which guarantees the emergence of a single-channel Kondo 
model at low energies. Local gates placed at the nanowire- 
superconductor contacts can tune the channel transparency. 
The nanowires are long enough so that their energy-level spac- 
ing is negligible, (b) All-metallic setup (side view). The 
atomic point contact fits a single conduction channel, whereas 
the wider NS junction may contain multiple channels. When 
L > (, Tk may be maximized by increasing the tunneling 
amplitude at the atomic point contact relative to that at the 
wide junction. 



Can the charge Kondo effect proposed in this pa- 
per be measured in experimentally realizable NSN de- 
vices? In a preceding subsection we have estimated 
charge Kondo temperatures of ~ 10 — 20mK for the 
single-channel Kondo model and reasonable values of A, 
r/<5 and E c . This temperature is close to the lowest tem- 
peratures achieved in current dilution fridge refrigera- 
tors (< 20mK). Moreover, Tk may in principle be made 
higher by designing more transparent junctions, tuning 
A/E c closer to one or choosing materials with larger val- 
ues of A. Thus while it will be difficult to experimentally 
access the unitary regime T <C Tk, measurements con- 
ducted in the weak coupling regime T > Tk may be 
within reach. 

Yet in spite of the extensive experimental work on 
mesoscopic NSN systems, there is no report of any CKE. 
Some trivial reasons for this include (i) the effective tem- 
perature of the device being too high, (ii) the gate voltage 
not being tuned with sufficient prccission to the charge 
degeneracy point, and (iii) purposely breaking the charge 
degeneracy by coupling the island to a superconducting 



reservoir. A more fundamental reason might be that the 
NS junctions fabricated in typical experiments lead to 
multichannel Kondo problems with concomitantly lower 

Next we discuss the experimental platform and re- 
quirements to observe CKE in NSN heterostructures. 
First and foremost we need A > E Cl which precludes 
ultrasmall superconducting islands due to their domi- 
nant energy^ Second, the superconducting island must 
be placed in the vicinity of a gate electrode. In order 
to measure CKE, the resolution of the gate voltage must 



be such that N a 



N, 



(0) 



where N g 



C g (V g V g {0) )/e < T K /E C , 
2M + 1 is the charge degeneracy point 
(for any integer M), and V g — eNg jC g is the cor- 
responding gate voltage. For typical values of the gate 
capacitances (C g ~ lOaF) and the Kondo temperature 
(Tk/E c ~ 10~ 3 ), it follows that the gate voltage resolu- 
tion needs to be better than 10/xV. This is achievable in 
current experimental devices. If N g — N g ' > Tk/E c , the 
Kondo RG equations (Eq. (fTTJ)) ) arc cutoff at energy scales 

- E c (N g - A^ 0) ) > T K and the charge Kondo effect will 
not be fully developed. However, even in this case there 
should be observable fingerprints of Kondo correlations 
(these are discussed in Section l*VT|) . 

Another desirable trait of the superconducting island 
is that it have no spatial symmetries and that its linear 
dimensions L be smaller than the superconducting coher- 
ence length £. This ensures having only one conduction 
channel in the dot at any given energy, which increases 
the likelihood of producing a single-channel Kondo model 
at low energies. Moreover, smaller island sizes mean a 
smaller ratio A/<5 and thus a larger value of T* , which en- 
hances the ultraviolet energy cutoff for the Kondo prob- 
lem. In principle, having a small-sized superconductor 
might conflict with the requirement that A > E c . In 
practice, there appears to be enough room in parameter 
space to ensure that A > E c and L < C, are satisfied 
simultaneously. For instance, the best developed ma- 
terial is aluminum, whose BCS coherence length in the 
ballistic limit is Co — hvp/^A) ~ 1.5/im. For a mean 
free path of I ~ 15nm, the actual coherence length is 
C = (HD/A) 1 / 2 ~ (Co0 1/2 - 150nm. Thus an Al island 
of L ~ 80nm connected to a normal metal with junc- 
tion area A ~ (80nm) 2 and e ~ 12 (dielectric constant 
of A10 x ) yields E c ~ 0.12meV < A ~ 0.2meV. Niobium 
is another attractive material, due to its higher critical 
temperature. Although its small BCS coherence length 
(Co ~ 40nm) requires relatively small and clean grains, 
the large dielectric constant of NbO x (e ~ 40) can pro- 
duce sufficiently small charging energies even for L ~ 
25nm and A ~ (25nm) 2 (E c ~ 0.4meV < A ~ 1.4meV). 
Tantalum and lead can also be interesting candidates, 
as they have BCS gaps that are larger than that of Al 
and coherence lengths that are longer than that of Nb. 
In Pb, Co — 90nm and C — 35nm for I ~ 15nm. With 
L ~ 30nm, A ~ (30nm) 2 and e ~ 25 (dielectric constant 
of PbO x ) we get E c ~ 0.4meV < A ~ lmcV. Note that 



these figures are merely representative. 

Last but not least, it is important that the NS in- 
terfaces have a small number of effective conduction 
channels. We have previously discussed how a single- 
channel Kondo model may be plausible even in wide 
metallic junctions with A 3> \ 2 F , where \p is the Fermi 
wavelength of the metal. In Fig. [3^ we sketch an al- 
ternative setup that combines a superconducting grain 
with semiconducting nanowires. This setup is moti- 
vated by recent success in manufacturing semiconduct- 
ing/superconducting nanocontacts i 39 ' 40 The large Fermi 
wavelength in semiconductors facilitates single-channel 
NS contacts without making the charging energy of the 
island larger than the superconducting gap. In addition, 
the possibility of tuning Schottky barriers in semiconduc- 
tor/metal interfaces by local gates creates prospects for 
measurable chargc-Kondo temperatures. 

Fig. [3)3 constitutes yet another all-metallic experimen- 
tal setup, which allows to reach sizeable Tk in NSN de- 
vices with L ^$> C (where electrons cannot coherently 
propagate from one junction to another). One of the 
junctions is made from an atomic point contact, while 
the other junction is wide enough to ensure that the 
charging energy of the device remains less than the su- 
perconducting gap. If the single channel at the point 
contact is much more strongly coupled to the supercon- 
ductor than each of the multiple channels at the wide 
junction, then the charge-Kondo temperature is close to 
that of the single-channel Kondo model. A theoretically 
straightforward way to achieve atomic point contacts is 
to use scanning tunneling microscopes (STM). STM tips 
operated in the point contact mode can^i capture An- 
dreev tunneling (which is altered by chargc-Kondo corre- 
lations), and recent STM experiments^ have succeeded 
measuring the density of states of ultrathin supercon- 
ducting Pb grains. An added challenge here is that the 
STM would have to operate at milliKelvin temperatures. 
Alongside milliKelvin STMs, mechanically controllable 
single-channel quantum point contacts are being devel- 
oped in superconducting systems 42 



D. Comparison to Other Charge-Kondo Effects 

Part of our motivation for exploring NSN systems is 
to establish a "quantum-dot counterpart" for the charge- 
Kondo effect that is believed to arise in compounds doped 
with negative-U molecules ("NU" systems). In this sec- 
tion we discuss similarities and differences, as well as ad- 
vantages and disadvantages of CKE in NU, NNN and 
NSN settings. Some of the highlights are summarized in 
Table El 

The ultraviolet (UV) cutoff A is the energy scale below 
which the Kondo model applies. It is lowest in NSN sys- 
tems, where T* <C A due to entropic arguments. This 
is partly why materials doped with negative-U centers 
reportedly show Tk — 10K, while we find Tk ~ lOmK 
in mesoscopic NSN heterostructures. Note that CKE in 



TABLE II: Comparison between different charge-Kondo effects (CKE) and the spin-Kondo effect (SKE). QD stands for semicon- 
ducting quantum dots, ui is the largest of a variety of extrinsic energy scales such as temperature, bias voltage, pseudo-magnetic 
fields and the frequency of an external perturbation. 





CKE in NSN 


CKE in NU 


CKE in NNN 


SKE in QD 


Two-level system 


N or N + 1 Cooper pairs 


N or N + 2 electrons 


N or N + 1 electrons 


f or 4- electron 


Ultraviolet cutoff (A) 


T*(5<T* <E C ) 


\U\ 


£ c (> S) 


<5(< E c ) 


Infrared cutoff 


UJ 


UJ 


max{w, 5} 


UJ 


fAx(A) 


- (r/s)f x 


~T/\u\ 


~ (r/*) 1/sl 


~r/E c 


"A,, (A) 


- (r/«J)/n 


= ^Ax(A) 





= v\±{A) 


Tk/S 


unrestricted 


~0 


needs to be large 


needs to be small 


Number of channels 


> 1 


1 


> 2 


> 1 



NNN systems does not suffer from entropic issues be- 
cause the ground state and the lowest excited state are 
both surrounded by similarly dense sets of single-particle 
states. 

The infrared (IR) cutoff u is the energy scale at which 
the renormalization of the Kondo couplings stops. In 
NSN and NU systems u> is the largest amongst the 
temperature T, bias voltage V, and pseudo-magnetic 
fields h z ,h x . The full Kondo effect is observable only 
if Tjc > u. In NNN systems, Kondo-like IR divergences 
originate from sums involving single-particle states in the 
dot. Consequently 8 too plays the role of an IR cutoff 
and the Kondo effect in NNN systems is observable only 
if Tk > max{<5, u>}. Tk > 5 is difficult to satisfy in semi- 
conducting dots. 

Concerning the Kondo couplings, in NU systems they 
are isotropic in pscudospin space and originate from a sin- 
gle energy level. In NSN and NNN systems, the Kondo 
couplings are anisotropic and sensitive to the multilevel 
energy spectrum of the dot. While the anisotropy be- 
comes gradually less pronounced at lower energies, it still 
plays an essential role in the quantification of Tk ■ More- 
over, multilevel effects can become interesting on their 
own (see Appendix). A unique feature of NSN systems 
is that the bare Kondo couplings are strongly depen- 
dent on A/E c , which can be tuned by magnetic fields. 
In SKE systems, the bare Kondo coupling is tipically 
smaller than in NSN systems because it contains a prod- 
uct of two small parameters: T/S and 5/E c . 

NNN systems can never have less than two Kondo 
channels, due to real spin degrees of freedom that are de- 
coupled from the charge pseudospin. In contrast, in NSN 
devices spin degrees of freedom are interwoven with the 
charge pseudospin (because the order parameter couples 
spin-up electrons with spin-down holes) and therefore the 
one-channel Kondo effect is in principle possible. Unlike 
NU and SKE systems, most NNN and NSN systems carry 
wide junctions with multichannel leads. However, multi- 
ple channels of a NSN system can be effectively mapped 
onto a one-channel Kondo model e.g. when the tunnel- 
ing matrix amplitude is independent of energy. Reaching 
a two-channel Kondo model in NNN systems with wide 



junctions is much less likely. 

In sum, CKE in NSN systems ought to be experimen- 
tally measurable much like ordinary SKE is observable 
in semiconducting quantum dots, provided that single- 
channel Kondo models can be engineered. 



VI. ZERO-BIAS CONDUCTANCE 

In previous sections we have described the charge- 
Kondo model that arises in mesoscopic NSN heterostru- 
cures at sufficiently low energies. In this section we in- 
vestigate its fingerprints in the low-temperature electri- 
cal transport. We focus on the elastic contribution to the 
zero-bias conductance, which is dominant at T < T* be- 
cause the superconducting gap freezes out inelastic quasi- 
particle excitations. Fingerprints of CKE can also get 
manifested in the low-temperature capacitance of super- 
conducting islands coupled to a single normal metallic 
lead^i^i 

The electrical conductance of mesoscopic NSN 
heterostructures has been widely discussed in the 
literature j 21 ' 29 ' 34 ' 44 " — nonetheless, there appears to 
be no report on many-body anomalies in the low- 
temperature elastic conductance. When bias voltages are 
weak, the current-voltage characteristics is encoded in 



I a =g aa 'V a {ae{L,R}), 



(25) 



where V a and I a = (I a ) = iesgn{a)([H c s, N a }} are the 
voltage drop and the current across the a-th junction, 
respectively. N a denotes the electron number operator in 
the a-th lead and ri c s is given by Eq. ([7]). Also, sgn(a) = 
1(-1) for a = L(R). 

The zero-bias conductance matrix in Eq. (|25[) can be 
expressed in terms of the Kubo formula: 



g aa , = lim — 







dTe lUT {T T i a {T)i a (0)), (26) 



where /3 = l/(fcsT) and T T is the time-ordering operator. 



10 



It is straightforward to obtain I a = I" in „ io + Ip air , where: 
4"„ g io = -iesgn(a) J^ X T' J2 S z c\^c a ^ c + h.c. 
7p air = iesgn(a)^A" Q J^ S + ac aka c a , k ^ a + h.c. 



a' k,k', cr 



(27) 



(^single) and (/pair) arc the single-particle and Cooper- 
pair tunneling contributions to the current, respectively. 
a 7^ a' in the expression for /single simply recognizes the 
fact that local (a = a') single-particle processes are nor- 
mal reflections that do not contribute to the current. We 
treat the leads as non-superconducting ((c a w\-c a k'i.) = 
0), non-magnetic ((c^ k ^c a k'^) = 0) and infinitely long 

(( c ak<r c "kv) oc (5k, k') conductors. With this in mind, the 
combination of Eqs. (f2"6]l and (|27|) yields 



Qaa' 






Gec + Gca Gec — Gca 
Gec — Gca G a + Gec + Gca 



where 



(28) 



G A {R) = 2G f dE(-d E f)2[^X L A L{RR \E)} 2 
Gca = G J dE(-d E f)2lv\i R (E)} 2 
Gec = Go I ^(-de/^A^tE)] 2 , (29) 



Go = ne 2 /h and f(E) = {cxp{f3E) + l)" 1 . G A {R} is the 
local Andreev conductance across the left (right) junc- 
tion, Gca is the non-local or crossed Andreev conduc- 
tance, and Gec is the elastic, single-particle co-tunneling 
conductance! 45 ' 47 The probability for an Andreev re- 



vO!a'\2 



flection process is proportional to (A" a ) + (A^ a ) 
2(A" Q ) 2 , whereas the probability for a co-tunneling event 
is ex (Aj^) 2 . The derivation of Eq. (g8|| is standard 4 ^ 
except for the following two details. First, A|| t ±(E) is 
energy-dependent via the Kondo RG flow of Eq. (|T9l . 
and diverges as E — > 0. Second, the current oper- 
ator in Eq. (|27[) contains S z and S ; we have used 
(%S z (t)S z (0)) = 1/4 and (T t S+(t)S~(0)) = 1/2 for 

T > 0. 

The matrix structure of g aa ' in Eq. (|28[) is simple to 
interpret. Regarding diagonal matrix elements, each lo- 
cal Andreev reflection contributes twice as much as each 
co-tunneling or crossed Andreev reflection event because 
the former involves two electrons crossing the same junc- 
tion while in the latter processes the two electrons cross 
different junctions. This fact is reflected by the extra 
factor of two in the first line of Eq. (f2"9"|. Concerning 
the off-diagonal matrix elements, electron co-tunneling 
and crossed Andreev reflection contribute with opposite 
signs because in the later case it is a hole (rather than a 
electron) that tunnels out of the superconducting island. 



Even though Eq. (|28|) formally agrees with expressions 
derived in previous studies of NSN heterostructuresj 4 ^— 
there are conceptual novelties when the superconducting 
island behaves as a pseudospin 1/2 quantum impurity. 
For instance, Gec/Gca — > 1/2 as the temperature is de- 
creased, irrespective of the strength of electron-electron 
interactions in the superconducting island. Therefore, 
crossed Andreev reflection and electron cotunneling do 
not cancel each other in the off-diagonal matrix elements 
of Eq. (|28[) . Other differences will become clear below. 

In steady state Ir = II = I and the electrical response 
to the total voltage drop Vr + Vl = V across the device 
can be evaluated by inverting Eq. (|28|) : 



r - dI _ r , ( g a + G A )G C A + G L A G\ 
G= dV- GEC+ G\ + G« + AGca ■ m 

Kondo correlations become apparent when we evaluate 
Eq. dSUJ) for T < T* . Let us suppose that the two NS 
junctions can be mapped onto a single- channel Kondo 
model. Combining Eqs. <(T3j) , flU]) and ((30]) it follows 
that 



G 



Gotrt 



L L R 



(«i + *l) 2 



dE(~d E f){(v\ ll ) 2 + 2(v\ ± ) 2 



(31) 
where X\\j_(E) obeys Eq. (IT9"1) . 

When A_l(T*) = A||(T*) (i.e. E c ~ 0.5A) and T* > 
T > T K , the integral in Eq. (J3TJ yields 



G — Go~ 



.2 .2 
l L l R 



m +& fog 2 (T/T K y 



(32) 



For T < Tk, the Kubo formula expression breaks 
down and the conductance converges to 2 - G = 
2(e 2 / h)4t 2 L t 2 R / (t 2 L + t 2 R ) 2 . The same results arise also 
in negative-U molecules with charge degeneracy,— in 
normal-metallic quantum dots with charge degeneracy 4 - 
and even in semiconducting quantum dots with spin 
degeneracy^ 

The conductance for X ± (T*) ^ A||(T*) at T* > T > 
Tk is displayed in in Fig. Q. G(T) increases with 
E c /A, as expected from Fig. ([T]). Since the anisotropy of 
the bare Kondo couplings disappears at sufficiently low 
temperatures, G(T <C Tk) is the same irrespective of 
Ax(T*)/A||(T*). 

The conductance curves of Fig. Q have been neither 
predicted nor measured in previous studies of NSN sys- 
tems. For instance, Rcf. |3J] neglected the inter-junction 
electron coherence (i.e. Gca = Gec — in Eq. (|3U)) ) 
and concluded that the zero-bias conductance through 
the superconducting grain will be independent of tem- 
perature, as long as T -C T* and the gate voltage is 
tuned to the charge degeneracy point. Contemporary 
experiment a 28 ' 29 appeared to agree with this theoretical 
prediction, although the temperature dependence of the 
zero bias peak was not analyzed in detail. Yet accord- 
ing to our theory, the zero-bias conductance could dis- 
play traces of a two-channel Kondo effect in islands with 
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FIG. 4: Low-temperature (elastic) zero-bias conductance 
(Eq. ([31}) across the NSN heterostructure, in the single- 
channel Kondo model. The conductance is normalized to 
its value at T = T*. Curves with red squares, black cir- 
cles and blue triangles have T K = 0.004T*, T K =± 0.01T* and 
Tk — 0.05T*, respectively. For simplicity we ignored g±,\\ 
in Eq. I|20[l and evaluated Tk analytically neglecting third 
order terms in Eq. JlQj) . We chose A = 10005. Inset: The 
solid line represents G(T)/G(T*) for E c = 0.5A, as calculated 
from Eq. (|31|) . Because E c = 0.5A yields the isotropic Kondo 
model (Ax(T*) = A||(T*)), the solid line matches well with 
Eq. (1321) (dashed line). 



L ^> Q . This is the case in the experimental setup of 
Fig. (3}d, which should host a low-temperature enhance- 
ment of the Andreev reflection at the point contact, ac- 
companied by a reduction of the Andreev reflection in 
the wide junction. 



VII. CHARGE-KONDO EFFECT IN THE 
STRONG TUNNELING REGIME 

Thus far all our calculations have concentrated on the 
weak tunneling regime, where T -C 5. In this section 
we turn to highly transparent NS junctions (T 3> 8) 
and determine the fate of CKE when charge in the su- 
perconducting island is strongly fluctuating. On one 
hand, in the weak tunneling regime the chargc-Kondo 
temperature is enhanced by more transparent normal- 
superconducting junctions (Section [Vj. On the other 
hand, when charge in the dot fluctuates strongly it is 
no longer appropriate to regard the superconductor as a 
two-level system. Which of these two conflicting trends 
prevails? It turns out that the answer to this ques- 
tion depends on electron-electron interactions in the lead, 
which we have ignored up until now. Strong tunneling 
and electron-electron interactions render bosonization as 
the most convenient method to employ in this section. 
In particular we combine and suitably modify the ap- 
proaches from Rcfs. [501 and [5T| . 



We consider a one-dimensional normal metallic wire 
connected to a superconductor (which need not be one- 
dimensional). The lead occupies the negative half- axis 
(x < 0), the NS interface is at x = 0, and the supercon- 
ductor is located at x > 0. We write the low-energy and 
long-wavelength Hamiltonian of this system as 



U = H 



lead 



Hc + Hb, 



(33) 



where 



"Hlcad = iVF J2 dx (*LA*£* ~ ^L^*^] 

_ J—oo 



is the Hamiltonian corresponding to the lead, 

J-co „ 

is the charging energy and 



(34) 



H 



Ra 



N n 



(35) 



n B = aV B J2 *£(0)M0) + aA B [* t (0)* i (0) + h.c] 

(36) 
is a boundary interaction (a is the lattice spacing). In 
Eqs. (f3"4l) . ([33]) and (|3l))) we have expanded the orig- 
inal fermion operators as 5 , cr (a;) ~ e~ lkFXl $>La{x) + 
e lkFX ^R cr {x) 1 where Uf is the Fermi wavevector and '4/ l<j 
(*& Ro) describes left-moving (right-moving) fcrmions. In 
Eq. (|3"B")) we have exploited total charge conservation in 
the dot-plus-lead system^ so as to rewrite the charging 
energy in terms of the charge density in the lead (in- 
stead of in the dot) . Eq. (|31))) follows from integrating out 
the superconductor at energy scales that are lower than 
the BCS gap A. This integration results in a boundary 
term^i which encodes normal scattering (ex Vb -C A) as 
well as Andreev scattering (ex Ag). Both Vb and As 
(which can be taken to be purely real in our case) de- 
pend on microscopic details of the superconductor such 
as its Fermi velocity, A and S. Also note that quasiparti- 
cle transmission into the superconductor is forbidden at 
energies below the gap. 

We bosonize Eq. (|3"3")) in the standard manner^ 



Vt --^ 5 [r$ c (a;)-e c (a;)+r(T* s (a;)-crG 3 (x)] 



* rCT (x) 



(37) 
where r\ a is a Klein factor (rj-fTJi = —fJiVti V<rVa — 1)> r = 
+ (— ) for right- (left-) moving fermions and a = +(— ) 
for spin-up (-down) fermions. $ c ($ s ) and 6 C (O s ) are 
the usual charge- (spin-) bosons. 
It follows that 



H 



lead 



VF 

2- 



E 



d.v 



K a (d x Q a ) 2 + -^(d x $ a ) 2 



(38) 

where K c (K s ) is the Luttingcr parameter in the charge 
(spin) sector (K s = 1 due to spin rotational symmetry). 
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Similarly, 



Hc=E c [V2$ c (0)/ir + N g 



(39) 



where ^/2^ c (0)/tt equals the number of electrons con- 
tained in x € [—oo,0] up to a constant (we take 
* c (-oo) = 0). 

In order to bosonizc Us, we use boundary condi- 
tions that arc associated with perfect Andrcev reflec- 
tion at the NS interfaced *Hf(°) = ~**U( ) and 
^Ri(0) = ^Lf(O)- These boundary conditions provide 
a reasonable approximation for highly transparent NS 
interfaces with small Fermi surface mismatches. Fur- 
thermore, perfect Andrecv reflection constitutes a renor- 
malization group fixed poinl^i when E c = Vb = 0. In 
bosonic language, the Andrcev boundary conditions can 
be translated as \/2Sc(0) = 7r/2 and v/2$«(0) = n. Con- 
sequently the boundary interaction at the Andreev fixed 
point reads 



H, 



-(2/7r)F B cos[\/2$ c (0)] 



(40) 



which describes single-particle backscattcring (normal re- 
flection) at the interface. Terms proportional to As do 
not appear in Eq. (|40|) . which is consistent with the fact 
that boundary Andreev reflection is never a relevant per- 
turbation at the Andreev fixed point ^ 

Next, we determine how E c and Vb evolve as the tem- 
perature is lowered. Assuming that their bare values are 
small compared to A, the leading order RG flow equa- 
tions read 



dE c 
dl 



= E n 



dV B 
dl 



= (1 - K C )V B 



(41) 



where I = ln(A/T). Thus E c is a relevant perturbation, 
and Vb is relevant, marginal or irrelevant depending on 
whether electron-electron interactions in the lead are re- 
pulsive (K c < 1), non-interacting (K c = 1) or attractive 
(K c > 1), respectively. This is markedly different from 
the non-superconducting system analyzed in Ref. [50J |. 
where dVsjdl = (1 — K c /2)Vb and backscattering is a 
relevant perturbation for non-interacting electrons. The 
underlying reason for this qualitative difference is that 
a perfectly transmitted channel at a NN interface has a 
conductance of 2e 2 /h, whereas a perfectly transmitted 
channel at a NS interface has a conductance of 4e 2 //i 
(due to perfect Andreev reflection). 

Since either E c or Vb is certain to grow under RG, 
the Andreev fixed point is unstable. The nature of the 
new fixed point depends crucially on whether E c or Vb 
reaches strong coupling first. If E c wins, the charge boson 
at the boundary gets pinned to a value dictated by the 
gate voltage, ^/2<& c (Q)/ir — > —N g , and the ground state is 
non-degenerate. Since the gate charge N g is a continuous 
variable, the discreteness of charge is unimportant in this 
case and the system is nowhere equivalent to a two-level 
system. 

If Vb prevails instead, then the charge boson at the 
interface becomes pinned through ^/2^ c (0)/tt — > 2n, for 



any integer n. The discreteness of charge becomes impor- 
tant through the discreteness of n, although there are an 
infinite number of degenerate ground states. This energy 
degeneracy between different values of n is lifted by Tic 
(Eq. ([39]) '). Choosing N g appropriately (e.g. N g = 1) the 
lowest energy state becomes two-fold degenerate (n = 
and n = —1), where the two states differ from one an- 
other by two electrons. The remaining minima of the 
cosine potential are at energies > E c higher, where E c 
is the renormalized charging energy. Therefore in the 
strong coupling limit of Vb the ground state of Eq. (j3"3"|) 
behaves as a two-level quantum system. Rare tunneling 
events between the n = and n = — 1 charge states (in- 
duced by the kinetic energy terms in Eq. (|3"5)) ) are tanta- 
mount to the rare spin-flip events of a magnetic impurity 
that is weakly coupled to a bath of itinerant fermions. 
Those tunneling events become prominent as energy is 
lowered further, and ultimately culminate in the charge- 
Kondo effect. 

In sum, for CKE to emerge in the strong tunneling 
regime it is necessary that Vb dominate over E c under 
renormalization group. In terms of the parameters of the 
original Hamiltonian, this condition reads 



> 



E c (0) \V B {0) 



K C <1, 



(42) 



where Vb(0) and E c (0) are the values of the backscatter- 
ing potential and the charging energy at energy scales of 
order A. In mesoscopic NSN devices Eq. (|42|) is difficult 
to satisfy unless electron-electron interactions in the lead 
are strongly repulsive and/or the bare charging energy 
is very small. The temperature at which Vb reaches the 
strong coupling regime is To = Ae~ lc , where l c is defined 
through V B (l c ) - A. Thus 



T ~ A(F b (0)/A)t^ <A. 



(43) 



If To ^> E c (l c ) = T c (0)e' c , inelastic charge excitations 
(i.e. transitions to n / 0, —1) are present at tem- 
peratures T) > T > E c (l c ). In this case, by anal- 
ogy with Ref. [50], the dimensionless Kondo coupling is 
v\ ~ E c (l c )/T and the Kondo temperature scales like 
T K ~ E c (l c ) exp(-l/i/A) « E c (lc) « T). If T < E c (l c ), 
inelastic excitations arc forbidden as soon as T < To and 
Tk could be a sizeable fraction of T). Yet because To is 
typically very small in weakly interacting leads (K c ~ 1), 
Tk is unlikely to be detectable. Prospects might be 
better in leads made from carbon nanotubes, where^ 
K c ~ 0.2 < 1 

In Section [V] we learned that the charge-Kondo tem- 
perature is unobservably low when T/S — > 0, though it 
increases rapidly as T/S grows. In this section we have 
seen that Tr- is very small when T/S 3> 1 as well, at least 
for Fermi liquid leads. Therefore, intermediate tunneling 
strengths are most suitable for the observability of CKE 
in NSN devices. 
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VIII. SUMMARY AND CONCLUSIONS 

Partly motivated by the recent experimental search 
for chargc-Kondo correlations in Tl-doped PbTe, we 
have developed a theory for the charge- Kondo effect 
in artificially fabricated, chaotic superconducting islands 
that are weakly connected to non-superconducting leads. 
We have focused on tunable superconducting chargc- 
qubits where the energy gap in the excitation spectrum 
is larger than the electrostatic charging energy. The 
low-energy transport and thermodynamic properties of 
these systems showcase Kondo-like behavior when the 
electrostatic energies of two charge states differing by 
a Cooper pair become degenerate. A single-channel 
anisotropic Kondo model is found to describe supercon- 
ducting grains of size smaller than the superconducting 
coherence length, which are coupled to normal conduc- 
tors either through atomic point contacts or semicon- 
ducting nanowircs. The same model could also apply 
for wider metal-superconducting contacts as long as the 
tunneling amplitude is independent of energy in a nar- 
row strip around the Fermi energy. For this model the 
Kondo temperature can reach ~ lOmK with aluminum 
grains and ~ 20mK with niobium grains, which is compa- 
rable to the temperature of the cryostat in current NSN 



devices. Higher Kondo temperatures might be acces- 
sible by engineering normal-superconductor contacts of 
intermediate transparency. Our proposal for the chargc- 
Kondo effect differs qualitatively from previous propos- 
als in non-superconducting single-electron devices and in 
valence-skipping compounds, and may be easier to ob- 
serve. 

Avenues for further theoretical research in NSN de- 
vices include understanding the transition from single- 
channel to multichannel Kondo models, examining the 
crossover from weak to strong dot-lead tunneling, char- 
acterizing the chargc-Kondo effect when the normal- 
superconductor junctions are disordered, and searching 
for a low-temperature many-body enhancement of the 
tunneling magnetoresistance when the leads are magne- 
tized. 



Acknowledgments 

I am indebted to I. Affleck for many helpful discus- 
sions. This research has been funded through a Junior 
Fellowship of the Canadian Institute for Advanced Re- 
search. 



1 For a review see e.g. D.L. Cox and A. Zawadowski, Adv. 
Phys. 47, 599 (1998). 

2 L.I. Glazman and M. Pustilnik, in Nanophysics: Coher- 
ence and Transport (eds. H. Bouchiat et al.), 427 (Elsevier, 
2005) 

3 K.A. Matveev, Zh. Eksp. Teor. Fiz. 99, 1598 (1991). 

4 L.I. Glazman and K.A. Matveev, Zh. Eksp. Teor. Fiz. 98, 
1834 (1990). 

5 K.A. Matveev, Phys. Rev. B 51, 1743 (1995). 

6 G. Zarand, GT. Zimanyi and F. Wilhelm, Phys. Rev. B 
62, 8137 (2000). 

7 A. Buxboim and A. Schiller, Phys. Rev. B 67, 165320 
(2003) 

8 K. Le Hur, P. Simon and L. Borda, Phys. Rev. B 69, 
045326 (2004) 

9 P. Kakashvili and H. Johannesson, Europhys. Lett. 79, 
47004 (2007). 

10 A. Taraphder and P. Coleman, Phys. Rev. Lett. 66, 2814 
(1991). 

11 P.S. Cornaglia, H. Ness and DR. Grempel, Phys. Rev. 
Lett. 93, 147201 (2004). 

12 M. Dzero and J. Schmalian, Phys. Rev. Lett. 94, 157003 
(2005). 

13 J. Koch, E. Sela, Y. Oreg and F. von Oppen, Phys. Rev. 
B 75, 195402 (2007). 

14 M.S. Laad, L. Craco and A. Taraphder, Europhys. Lett. 
88, 20008 (2009). 

15 D. Goldhaber-Gordon, H. Shtrikman, D. Mahalu, D. 
Abusch-Magder, U. Meirav and M.A. Kastner, Nature 
391, 156 (1998); S.M. Cronenwett, T.H. Oosterkamp and 
L.P. Kouwenhoven, Science 281, 540 (1998); R.M. Po- 
tok, I.D. Rau, H. Shtrikman, Y. Oreg and D. Goldhaber- 



Gordon, Nature 446, 167 (2007). 

16 T. Gramespacher and K.A. Matveev, Phys. Rev. Lett. 85, 
4582 (2000). 

17 E. Lebanon, A. Schiller and F.B. Anders, Phys. Rev. B 68, 
155301 (2003). 

18 D. Berman, N.B. Zhitenev, R.C. Ashoori and M. Shayegan, 
Phys. Rev. Lett. 82, 161 (1999). 

19 K.W. Lehnert, B.A. Turek, K. Bladh, L.F. Spietz, D. Gun- 
narson, P. Delsing and R.J. Schoelkopf, Phys. Rev. Lett. 
91, 106801 (2003). 

20 Y. Matsushita, H. Bluhm, T.H. Geballe ad IR. Fisher, 
Phys. Rev. Lett. 94, 157002 (2005); M. Matusiak, E.M. 
Tunnicliffe, JR. Cooper, Y. Matsushita and IR. Fisher, 
Phys. Rev. B 80, 220403(R) (2009). 

21 A.D. Zaikin, Physica B 203, 255 (1994). 

22 Y. Dubi, Phys. Rev. B 75, 094510 (2007). 

23 M. Tinkham, Introduction to Superconductivity, 2nd ed. 
(Dover, New York, 1996). 

24 J. von Delft and D.C Ralph, Phys. Rep. 345, 61 (2001). 

25 I.L. Aleiner, P.W. Brouwer and L.I. Glazman, Phys. Rep. 
358, 309 (2002). 

26 D.V. Averin and Yu. V. Nazarov, Phys. Rev. Lett. 69, 
1993 (1992). 

27 P. Lafarge, P. Joyez, D. Esteve, C Urbina and M.H. De- 
voret, Phys. Rev. Lett. 70, 994 (1993). 

28 T.M Eiles, J.M. Martinis and M.H. Devoret, Phys. Rev. 
Lett. 70, 1862 (1993). 

29 J.M. Hergenrother, M.T. Tuominen and M. Tinkham, 
Phys. Rev. Lett. 72, 1742 (1994). 

30 Y. Nakamura, Yu. A. Pashkin and J.S. Tsai, Nature 398, 
786 (1999). 

31 M. Houzet, D.A. Pesin, A.V. Andreev and L.I. Glazman, 



14 



Phys. Rev. B 72, 104507 (2005). 

G.A. Fiete, G. Zarand, B.I. Halperin and Y. Oreg, Phys. 

Rev. B 66, 024431 (2002). 

L.I. Glazman and M.E. Raikh, Pis'ma Zh. Eksp. Teor. Fiz. 

47, 378 (1988). 

F.W.J. Hekking, L.I. Glazman, K.A. Matveev and R.I. 

Shekhter, Phys. Rev. Lett. 70, 4138 (1993). 

A.C. Hewson, The Kondo Problem to Heavy Fermions 

(Cambridge University Press, Cambridge, 1993). 

P. Coleman, Introduction to Many-Body Physic s, 

http://www.physics.rutgers.edu/~coleman/mbody.html 

I. Garate and I. Affleck, Phys. Rev. Lett. 106, 156803 

(2011). 

A.O. Gogolin, A.A. Nersesyan and A.M. Tsvelik, 

Bosonization and Strongly Correlated Systems (Cambridge 

University Press, United Kingdom, 1998). 

J. Xiang, A. Vidan, M. Tinkham, R.M. Westervelt and 

CM. Lieber, Nature Nanotechnology 1, 208 (2006). 

Y.-J. Doh, S. De Francheschi, E. Bakers and L.P. Kouwen- 

hoven, Nano Lett. 8, 4098 (2008). 

J.Y.T. Wei, N.-C. Yeh, D.F. Garrigus and M. Strasik, 

Phys. Rev. Lett. 81, 2542 (1998) 

T. Nishio, T. An, A. Nomura, K. Miyachi, T. Eguchi, H. 

Sakata, S. Lin, N. Hayashi, N. Nakai, M. Machida and Y. 

Hasegawa, Phys. Rev. Lett. 101, 167001 (2008). 

E. Scheer, P. Joyez, D. Esteve, C. Urbina and M.H. De- 



voret, Phys. Rev. Lett. 78, 3535 (1997); C. Brun, I.-P. 
Hong, F. Patthey and W.-D. Schneider, [arXiv:1006.0"333V l 
(2010); M. Zgirski, L. Bre theau, Q. Le Mas ne, H. Pothier, 
D. Esteve and C. Urbina. larXiv:1105. 50391 (2011). 
M. Tinkham, J.M. Hergenrother and J.G. Lu, Phys. Rev. 
B 51, 12649 (1995). 

G. Falci, D. Feinberg and F.W.J. Hekking, Europhys. Lett. 
54, 255 (2001). 

G. Bignon, M. Houzet, F. Pistolesi and F.W.J. Hekking, 
Europhys. Lett. 67, 110 (2004). 

D.S. Golubev and A.D. Zaikin, Phys. Rev. B 76, 184510 
(2007). 

A. Levy Yeyati, F.S. Bergeret, A. Martin-Rodero and T.M. 
Klapwijk, Nature Physics 3, 455 (2007). 
G.D. Mahan, Many-Particle Physics, Physics of Solids and 
Liquids Series, 3rd ed. (Kluwer Academic/Plenum Publish- 
ers, New York, 2000). 

L.I. Glazman, F.W.J. Hekking and A.I. Larkin, Phys. Rev. 
Lett. 83, 1830 (1999). 

I. Affleck, J.-S. Caux and A.M. Zagoskin, Phys. Rev. B 62, 
1433 (2000). 

T. Giamarchi, Quantum Physics in One Dimension (Ox- 
ford University Press, New York, 2003). 
C. Kane, L. Balents and M.P.A. Fisher, Phys. Rev. Lett. 
79, 5086 (1997). 



J 



Appendix A: "First-Principles" Derivation of CKE in NSN Devices 



In Section ITVl we introduced a low energy effective Kondo model, Eq. ([7]), from which we inferred the Kondo 
effect by simply invoking standard renormalization group equations. While we provided physical arguments to justify 
the emergence of such a model at temperatures lower than T* , we did not present a rigorous derivation. The main 
objective of this Appendix is to validate Eq. ([7]) by deriving hallmark Kondo-like divergences from the full microscopic 
(or "first-principles") Hamiltonian, Eq. (Q]). 

The procedure we follow has been introduced elsewhere^ The task undertaken here consists of evaluating the 
amplitude for Andrccv reflection to fourth order in tunneling. Andrccv reflection transforms an itinerant electron into 
an itinerant hole, while adding a Cooper pair to the superconducting island. This can be interpreted as a process 
that flips the pseudospins of the superconducting island as well as of the incident fermion. The initial and final states 
for this process are 



|i)=4 t |0)|2M) ;|/)=c fc| |0)|2M + 2>, 



(Al) 



which arc eigenstatcs of Eq. (UJ in absence of tunneling. Tunneling induces \i) — > \f) transitions, and the transition 
amplitude Ai-^f may be computed from certain matrix elements of an effective Hamiltonian derived perturbatively 



from Eq. ^)M One arrives at Ai_»/ — A 



(2) 



A 



(4) 



>/ 



., where 



.4 



(2) 



.4 



(4) 



E 



(f\n T \n)(n\n T \i) 



E 

n 1: n 2 ,n 3 



Ei — E n 

(f\nT\n 1 )(n 1 \rlT\n2)(n 2 \n T \n 3 )(n 3 \'HT\i) 
(Ei — E ni )(Ei — E, l2 )(Ei — E n3 ) 



v {f\H T \n){n\H T \i) (2) 
2 4- [E. h -E n Y 



E 



\(f\n T \n)f 
(Ei - E n y 



and \n), \ni), ... denote intermediate states of energies E n , E ni ...in absence of tunneling. In addition, 



(2 



E 



E n — Ei 



= * 2 EE°fe) 



E c -i q - v/e^+A^ 3£ c + £ 9 + v^+A^ 



(A2) 
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where in the second equality we have assumed particle-hole symmetry and energy-independent tunneling amplitudes. 
Q(x) is the step function. 

vA^J.1 is identical to (T/8)f± of Section llVl and was first computed in Rcf. (34J. It is finite and non-singular for 



iW 



E c < A. Yet A\_^f, which has not been previously computed, hosts infrared (IR) logarithmic divergences for any 
E c /A. These IR divergences signal the onset of Kondo correlations. 

Evaluating A^ , is straightforward in principle but cumbersome in practice. For convenience we attach Tables 

IIIIIHVI which collect all possible intermediate state configurations for the first term in the expression of A^ t ■ Many 
of the individual amplitudes entering Tables IIIIIIIVI arc ultraviolet (UV) divergent in the infinite-bandwidth limit. 
Remarkably, all UV divergences in A^, end up cancelling one another after summing over all the amplitudes. This 
indicates that details of the energy spectrum at high energies are not important for the Kondo effect, and endows 
universality to results derived in this paper. In fact, the bulk contribution to j4i->/ originates from states with 

M < a. 

Following a numerical sum of all the amplitudes, we arrive at 



vA, 



-/ ± + 2— / ± /|,ln- 
o z " u 



r 2 



(A3) 



where g± is a dimcnsionlcss function of E c /A. To the present order approximation, Eq. (|A3|) can be recasted as 



vA^j ~ i/\± + 2(i/Aj.)(z/Aii)ln 



T * 



(A4) 



where 



^A± = -/i 
o 



£51 

Sfi 



r a 



and 



^ A n = t/ii 



TgnTfi A 



(A5) 



are the values of the Kondo couplings at the energy scale T* (see Eq. (|20|)) and we have neglected 0(T 3 ) terms. 
Eq. (|A4|) provides a first-principles proof for the second line of Eq. (fT!?|) and establishes the anisotropic Kondo model 
as an appropriate low energy theory for Eq. ([I]). Admittedly, the energy scale T* does not appear naturally in our 
zero-temperature calculation. Rather we introduce it in hindsight (i.e. with the knowledge that only at energies below 
T* can the superconducting island behave as a two-level system) and arrange the rest of the terms accordingly. The 
logarithmic terms in Eq. (|A5p indicate that the Kondo couplings begin to renormalize starting at energy scales of order 
A. g± is plotted in Fig. ([1]), where it is shown that 4th order tunneling events can cither enhance (if g± > 0) or deplete 
(if g± < 0) the Andreev reflection amplitude, depending on E c /A. We have not evaluated g\\, which would require 
additional Tables of intermediate state configurations. We expect g\\ ~ g± for E c ~ 0.5 A. Our pcrturbative analysis 
becomes insufficient when (r/(5)|gjj ~ 1, which occurs at E c ~ 0.8A for T/6 ~ 0.1. Of course, the breakdown of the 
perturbative regime does not preclude a Kondo effect; it just means that we cannot compute the Kondo temperature 
reliably using the weak-tunneling approach. 
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TABLE III: Virtual elastic processes to fourth order in single-particle tunneling, n, m and q are dummy variables to be summed 
over. The amplitudes labelled with a number are either UV and/or IR divergent. The amplitudes labelled with a letter are 
non-divergent. The amplitudes grouped with the same number or letter are identical to one another in presence of particle-hole 
symmetry. For each of these groups, half of the configurations involve electrons (cj) and half involve holes (c q ). For simplicity 
we have assumed the tunneling amplitude to be simply a constant; however, our results can be readily generalized to include 
more complicated situations. In the derivation of Fig. |T|) we use ^ —> v J f d£ and J^ n — > (1/8) f dn. N.B.: in these Tables 

(and in these Tables only) E n stands for V e n + A 2 . 
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c qrr \0hl t jl a \2M + 2} 


4t c ^l0)7^|2M) 


r, u n v n u' ! m &(-( -q ) 
Z (E C -E n )(i q -E n -E m )(i q + E c -E m ) 


5 


|0> 7 t f |2M + 2) 


4 t |0>7i t 7^|2M + 2> 


4 t 4j0)7i t i 2M ) 


u m v m vl&(( g ) 


(E a -E n )(E n +E m ){(, q +3E a + E n ) 


6 


|0) 7 t f |2M + 2) 


4 t |0)7l t 7^|2M + 2) 


4t4»7™j2A/) 


u n v„v^e(i q ) 


(E a -E n )(E n +E m )(S, q +3E c + E m ) 


■>■) 


|0) 7 t f |2M + 2) 


4»7i t 7Lr|2M + 2) 


4i4 t l0)7i t |2M) 


«m%<e({,) 




(.E c -E„)(.( q +E n +E m )(£ q +3E C +E„) 


10 


|0hl t |2M + 2) 


c 9 t|0)7i t 7i t |2M + 2) 


c q t4 t l0)7l t |2M + 2) 


u>nV n u^ n <r){ — £ q ) 


(E a -E n )(£ q -E n -E m )(E a +£ q -E n ) 


2 


|0}7i t |2M + 2) 


c,tl0>7it7ltl 2M + 2 ) 


|0)7^ t l 2M ) 


U n V n U~ n O(—t; q ) 


(E c -E n )(( q -E n -E m )(E c -E m ) 



E cvj.cJa|0>7i,- ff |2M + 2) 



4 CT |0)7lt7l- CT |2M + 2) 



|0} 7 t |2M + 2) 



v. m y m v„e(t q ) 



(E c -( q -E n )(-i q -E n -E m )(E -E m ) 



8 Cfc ^4» 7 : t |2A/ + 2> c fe ^|0) 7 ^ 7 t t |2A/ + 2) 



|0} 7 l t |2A/ + 2) 



u n v n u r 



,&(£„) 



(E c -£ q -E Tl )(-E m -E rh )(E -E n ) 



B c fc ^4 CT |0) 7 t j _ CT |2M + 2) 4t c fc^4 CT |0)7i,- CT 7lj2M + 2) 4 fCfc ^|0) 7 ^|2M + 2) 



^ m f m <e(eq) 



'(^-{,-B n )(-5,-E„-E m )(B c -B m ) 



9 


c & u4tl0>7lj2M + 2) 


owi\0)lUlU\2M + 2) 


C^4 t |0)7^|2M + 2) 


UnV-nU^&^q) 


(E c -!- q -E n )(-E n -E m )(E c -E„) 


10 


c^4j0) 7 l t |2M + 2) 


4j0)7l t 7^|2M + 2) 


|0} 7 l t |2M + 2) 


u n v n v'^e(s q ) 


(E a -H q -E n )(-(, q ~E n -E m )[E c -E n ) 


11 


c*';4tl0)7™j2M + 2) 


4 t |0)|2M) 


|0> 7 l t |2Af + 2) 


U n V n U m S(t; q ) 


(E a -Z q -E n ){-i q )(E B -E m ) 


H 


^4»7ltl 2M + 2 > 


4t c ^4>>i 2M > 


4t c ^l0)7lj2M + 2> 


U n D n U m 9(^q) 




(E -S q -E„)(-£ q )(E c -E m ) 


12 


c*u4tl0)7lj2M + 2> 


4 t |0)|2M) 


4t4tl0>7ij2M) 


v.v.v-nv'^&^g) 


(E c -i q -E n )(-Z q )(-Z q -3E c -E m ) 


13 


c fc ^4 CT |0) 7 t,_ a |2M + 2) 


4.|0)7^ t 7i,-j2M + 2) 


4<,4 t l0)7i,- CT |2M) 


(Bc-€<j-Br 1 )(-?,--B re --B m )(-€,-3Bc--Br l ) 


C 


^4J )7ltl 2M + 2) 


4j0)7i t 7i t |2M + 2> 


44 t l0)7l t |2M) 


«»««<8{S,) 


(B c -5,-E„)(-5,-i5 n -i5 m )(-5 g -3Ee--B m ) 


14 


c^4J )7ltl 2M + 2> 


4 t c fc ^c^| )|2Af) 


4 t 4j0)7lti 2M ) 


« n «n«4e(e„) 


(i5 c _ 5g _ B „)(_^)(_ 5g _ 3Bc _ Bm ) 


15 


Cfc';4tl0)7™j2M + 2) 


4 t c^4 t |0> 7 t i7 j r j2M + 2) 


4t c fe'il0)7lj2Af + 2} 


n„i.„i.^„e(e,) 


(Bo-€<j--B re )(?, + -B„+B m )(Bc--E„) 


13 


c^4j ) 7 t t |2M + 2) 


4t c feu4j0)7it7lj2M + 2> 


4 t 4ii0)7i t i 2M ) 


« m « m «iie(5,) 


(E c -( q -E n )(-i q -E n -E m )(-H q -3E c -E n ) 


13 


c&u4tl0)7ij2M + 2> 


4 t c fe a4tl0)7ni7^|2M + 2> 


4t4tl0)7nj2M) 


Um.v m v^e(Sq) 


(Ec-£ q -B n )(-$ q -E n -E m )(.-£ q -3E c -E n ) 


C 


c*u4tl0)7lj2M + 2) 


4 t c fc ^4 t |0) 7 ^ 7 lj2M + 2) 


ct t 4tl0>7^|2M) 


« n %^e(5,) 


(i3 c -4,-B„)(-e,-i3„-B m )(-e,-3B c -B m ) 


7 


c fc ^4j0)7l t l 2M + 2) 


c fc 'i|0>7l t 7lj2M + 2) 


Cfe'i4 t |0>7lj2M) 


u T1 i; ri ii TTi 0(^q) 


(B c -5 g -B„)(-B„-B m )(B c -B m ) 


7 


^4tl )7lj2M + 2) 


C^|0>7^7i t |2M + 2) 


|0)7l t |2M) 


U^ll^li^Gf^g) 


(Bc-€,--E„)(-E„-B m )(i3 c -i5 m ) 



17 



TABLE IV: Continuation of Tabled] 



Labc 


M 




M 


l»s> 


T4l ~A 

Contribution to A\_^ , (in units of t ) 


16 


c fe 'i4>hlj2A/ + 2) 


Cfe'j 


4 t 4 CT |0)7l,- CT 7lj2M + 2) 


Cfc^4tl0)7^|2M + 2) 




16 


c^4 t |0)7^|2M + 2) 


Cfc 


^4t c ^!^7l CT 7i4.|2A/ + 2) 


c^4 t l0>7^|2M + 2) 


A (E a -E n )^(i q -E n -E m ) 


7 


c^4tl^7^|2M + 2) 




4tl0)7lWntl 2M + 2 ) 


4t c 9 t|0)7i t |2M + 2) 


u m v m vf l e{-^ q ) 


(B c -B„)(B„+B m )(5,+B c -B m ) 


17 


c^4 t |0)7lj2M + 2> 


Cfe 


'^4t4tl^7^7lj2M + 2) 


c^4 t l0)7^|2M + 2> 


« m i' m i'^e(e,) 


(E c -E n )(( q +E n +E m )(E c -E m ) 


18 


c^4 t |0>7^|2M + 2) 




c fe u4t c 9t!0)!2M + 2) 


c^4 t |0)7^|2M + 2) 


» m « m ^e(-y 


(B»-B»)£ a (B e -Bm) 


H 


^4>hlj2M + 2) 




Cfe'i4t c 9tl0)|2M + 2) 


4t c 9 tl0hLr|2M + 2) 


»m%<8(-{,) 




(-Bc--B„)«,(i5 c +5 g -B m ) 


9 


c/bu4 t l0)7^|2M + 2) 




4tl0)7l i 7 7 1 ; it |2M + 2) 


4 t ^|0>7lj2M + 2) 


»„»„»^e(-£,) 


(B c -B„)(B„+B m )(B c +4,j-B„) 


18 


c/bu4 t |0)7^|2M + 2) 




c^4 t 4>>i 2M > 


c^4 t |0)7lj2M + 2) 


Un%M m B(^q} 


(S c -S„)£«(-Ec—E m ) 


19 


c*/ + cl t |0>7^|2M + 2> 




^^4t4J )i 2M ) 


4 t 4»%Ui 2M > 


Kjl1'Tl''^(^) 


(E c -B„)e,(e g +3B C + B m ) 


20 


c^4 t |0>7^|2M + 2) 


Cfe 


>AA»^U\ 2M + 2 ) 


4 t 4tl0)7nJ 2 ^) 


»m» m «ie(f,) 


(Bc--Bn)(5g + Er 1 +-B m )(5,+3Ec+E re ) 


E 


c^4 t |0>7^|2M + 2) 


Cfc 


^4t c 9 -l0)7^7ij2M + 2) 


c[, t c 9CT |0)7 t I[T |2M} 


^ (E a -E n )(e, q -E m -E n ){E c +Z q -E m ) 


5 


c*/ 4 ct t |0>7^|2M + 2> 




4tl0)7V:ntl 2M + 2 ) 


4t4tl0)7lj 2 ^) 


u. m v m v^e(S q ) 


(E c -E n )(E n + E m )(i q +3E c + E n ) 


D 


c^4tl0>7^|2M + 2) 


Cfc'j 


4 t 4»7l,- CT 7, 1 j2M + 2) 


4t4.l0)7l,-j2M} 


r, UnVKV^eUq) 

Z (E c -E n )(H q + E n + E m )(H q +3E c +E m ) 


17 


c/bu4 t |0>7^|2M + 2) 


Cfe 


^4 t ^|0>7^7l 4 .|2M + 2) 


c & u4 t l0)7lj2M) 


u n v n u^ n S( — £ q ) 


(E c -E n )(i q -E n -E m )(E c -E m ) 


15 


c^4tN)7^|2M + 2) 


Cfc 


<+4 t ^|0h;U7, t j2M + 2> 


4t c ^l0)7lj2M + 2) 


u n v n u rn S( — ^ q ) 


(E c -E n )tt q -E n -E m )(E a +t; q -E n ) 


21 


c^4tl0>7^|2M + 2) 




4tl0)7^7ltl 2M + 2 ) 


4t4>>7Ul 2M > 


u„v n v^e(i q ) 


(B c -B„)(-E„-B m )(-3Ec-?,-E m ) 


D 


c fc '4.c, CT |0)7L|2M + 2) 




c 9CT |0>7L7l t l 2M + 2 ) 


|0) 7 l t |2A/ + 2) 


(Zq-3E c -E n )(( q -E n -E m )(E c -E m ) 


A 


c fc '4,c 9CT |0} 7 L|2M + 2) 


Cfc 


j.c gCT 4 f |0>7L7Lj2M + 2) 


c^4 t l0)7^|2M + 2) 


(e<j-3Be--B„)(«,-E„--B m )(Bc--E m ) 


3 


c^c 9i |0) 7 tj2M + 2) 




c,4.|0)|2M + 2) 


|0) 7 ^ t |2A/ + 2) 


^m^m^^^C - £q) 


(£q-3E c -E„)£ q (E c -E m ) 


14 


c^c,;|0) 7 ^|2M + 2) 




c ql \0)\2M + 2) 


4t c ^|0)7lj2M + 2) 


U m ^7Ti'U n t ^ , ( — ^ q ) 


(e, q --AE € -E n )i q ^ q + E c -E m ) 


19 


c^c gt |0) 7 ^ t |2M + 2> 




Ck>iC qt cl t \0)\2M + 2) 


W=!tl0>7^|2M + 2> 


'U 7 nV rn U n @( — £ tq ) 


(£ q -3E c -E n )£ g (E c -E m ) 


12 


Cfc'iC« t |0) 7 i t |2M + 2) 




cyic qt c{ t \0)\2M + 2) 


c 9 t4 t l0)7^l 2M + 2 ) 


UmVm.U n O( — £ q ) 


(€ q -3E e -E n )( q (i q +E e -E m ) 


13 


c fc '4.c, CT |0) 7 L|2M + 2) 




c q „\0}~fLlii\2M + 2} 


4t c 9-|0>7L|2M) 


O UmV m U^B(-i q ) 

^l,i q -3E c -E n )(S, q -E n -E m )(i q +E c -E n ) 


13 


c fe /;C, CT |0) 7 L|2M + 2) 


Cfc 


lC q ,cl r \0) 1 L 1 l l \2M + 2) 


c 9CT 4 t |0>7L|2M) 


r> Um»m"i8(-{fl) 
Z («,j-3Bc~E, 1 )(e,--B„-i3 m )(«<j+£;c-i3n) 


5 


Cfc'j.c 84 |0>7^|2M + 2> 




owi\0hl t lU\2M + 2) 


Cfea4tl0)7lj2Af) 


W m f m W ^ © ( — £ q ) 


^ q -3E c -E n )(-E n -E m )(E c -E n ) 


20 


c^c 9t |0) 7 t t |2M + 2) 




c 9 tl0>7l t 7l lt |2Af + 2) 


|0>7l t |2M) 


UmVmUn@( — £ q ) 


(£ q -3E c -E n )(t Lq -E n -E m )(E c -E n ) 


■•>■) 


c^c, 4 |0) 7 ^|2M + 2> 


Cfe 


';^4tl0>7^7^|2M + 2) 


c fe ';4 t |0)7lj2Af) 


'^■m'"m'^ n ®{ — £,q) 




(€ q -3E c -E n )(£ q -E n -Em)(Ec-En) 


5 


C*'4.Cat|0>7»t|2M + 2> 




c^|0>7^7lt|2M + 2) 


|0)7l t |2M) 


U m V m V- n © ( — £, q ) 


^ q -3E c -E n )(-E n -E m )(E c -E n ) 


6 


c^c gt |0) 7 + f |2M + 2) 




<Vi|0hVitl 2M + 2 > 


c^4 t |0)7^|2M + 2) 


u m'^m'^n®{ — £,q) 


(i q -3E c -E rl ){-E n -E m )(E c -E m ) 


21 


c^c, 4 |0) 7 ^|2M + 2> 




c fc ^|0)7Vlj 2M + 2 ) 


|0} 7 ^ t |2A/ + 2) 


UmVm'U' n O( — ti q ) 




(£, q -ZE c -E n )(-E n -E m )(E c -E m ) 


C 


C*^Cat|0>7»t|2M + 2} 




c 9 tl0)7l t 7l t ! 2M + 2 ) 


4t c 9 T|0>7^tl 2M + 2 > 


UmVm.U n O( — £ q ) 


(.t, q -3E c -E„)(,$ q -E n -E m )(( q +E c -E m ) 


C 


c^c 9i |0) 7 tj2M + 2) 


C/b 


'i^4 t |0)7^7^|2M + 2) 


4t c ^l0)7ij2M + 2) 


Urn Vjn U^ @ ( — £ q ) 


{Z q -3E c -E n )(£ q -E n -E m )(£ q + E c -E 1Tl ) 



